GENERALIZED FRACTIONAL CALCULUS WITH 
APPLICATIONS TO THE CALCULUS OF VARIATIONS 

TATIANA ODZIJEWICZ, AGNIESZKA B. MALINOWSKA, AND DELFIM F. M. TORRES 

Abstr act. We study operators that are generalizations of the classical Riemann- 
Liouville fractional integral, and of the Riemann— Liouville and Caputo frac- 
tional derivatives. A useful formula relating the generalized fractional deriva- 
tives is proved, as well as three relations of fractional integration by parts that 
change the parameter set of the given operator into its dual. Such results 
are explored in the context of dynamic optimization, by considering problems 
of the calculus of variations with general fractional operators. Necessary opti- 
mality conditions of Euler-Lagrange type and natural boundary conditions for 
unconstrained and constrained problems are investigated. Interesting results 
are obtained even in the particular case when the generalized operators are 
reduced to be the standard fractional derivatives in the sense of Riemann- 
Liouville or Caputo. As an application we provide a class of variational prob- 
lems with an arbitrary kernel that give answer to the important coherence 
embedding problem. Illustrative optimization problems are considered. 



1. Introduction 

Fractional calculus studies derivatives (and integrals) of non-integer order. It is 
a classical mathematical field as old as calculus itself [25]. During almost 300 years, 
fractional calculus was considered as pure mathematics, with nearly no applica- 
tions. In recent years, however, the situation changed dramatically, with fractional 
calculus becoming an interesting and useful topic among engineers and applied 
scientists, and an excellent tool for description of memory and heredity effects |30| . 

One of the earliest applications of fractional calculus was to construct a com- 
plete mechanical description of nonconservative systems, including Lagrangian and 
Hamiltonian mechanics [42l [43] . Riewe's results {42J [43] mark the beginning of the 
fractional calculus of variations and are of upmost importance: nonconservative 
and dissipative processes are widespread in the physical world. Fractional calculus 
provide the necessary tools to apply variational principles to systems character- 
ized by friction or other dissipative forces, being even possible to deduce fractional 
conservation laws along the nonconservative extremals [18] . 

The theory of the calculus of variations with fractional derivatives is nowadays 
under strong current development, and the literature is already vast. We do not 
try to make here a review. Roughly speaking, available results in the literature 
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use different notions of fractional derivatives, in the sense of Riemann-Liouville 
QH [37], Caputo [201 [321 [36], Riesz [U Ell E], combined fractional derivatives 
[2"71 12"81 [33], or modified/generalized versions of the classical fractional operators 
[3j [T3l HH EU [38], in order to describe different variational principles. Here we 
develop a more general perspective to the subject, by considering three fractional 
operators that depend on a general kernel. By choosing special cases for the kernel, 
one obtains the standard fractional operators and previous results in the literature. 
More important, the general approach here considered brings new insights and give 
answers to some important questions. 

The text is organized as follows. In Section[2]the generalized fractional operators 
Kp, Ap and Bp are introduced and basic results given. The main contributions of 
the paper appear in Section[3j we prove a useful relation between Ap and Bp (The- 
orem EHJ, several formulas of integration by parts that change the parameter set P 
into its dual P* (Theorems 13.41 137H and I3T5)) . and new fractional necessary optimal- 
ly conditions for generalized variational problems with mixed integer and fractional 
order derivatives and integrals (Theorems 13. 1 H f3 . 171 and l3.22j) . We see that even for 
an optimization problem that does not depend on generalized Riemann-Liouville 
fractional derivatives, such derivatives appear naturally in the necessary optimality 
conditions. This is connected with duality of operators in the formulas of inte- 
gration by parts and explains no-coherence of the fractional embedding [T3. . This 
is addressed in Section 13.41 where we give an answer to the important question 
of coherence, by providing a class of fractional variational problems that does not 
depend on the kernel, for which the embedded Euler-Lagrange equation coincides 
with the one obtained by the least action principle (Theorem I3.26p . Finally, some 
concrete examples of optimization problems are discussed in Section [U 



2. Basic notions 

Throughout the text, a denotes a positive real number between zero and one, 
and diF the partial derivative of a function F with respect to its ith argument. 
We consider the generalized fractional operators Kp, Ap and Bp as denoted in [2]. 
The study of generalized fractional operators and their applications has a long and 
rich history. We refer the reader to the book [26] , 

Definition 2.1 (Generalized fractional integral). The operator Kp is given by 

t b 
Kpf(t)=P J k a (t,r)f(r)dT + q J k a (r, t)/(r)rfr, 

a t 

where P — (a,t,b,p,q) is the parameter set (p-set for brevity), t g [a, b], p,q are 
real numbers, and k a (t,r) is a kernel which may depend on a. The operator Kp 
is referred as the operator K (K-op for simplicity) of order a and p-set P. 

Theorem 2.2. Let a G (0,1) and P = (a,t,b,p,q). Ifk a {-,-) is a square-integrable 
function on A = [a, 6] x [a, b], then Kp : Li ([a, 6]) — > L2 ([a, b]) is a well defined 
bounded linear operator. 



Proof. Let a 6 (0, 1) and P = (a, t, b,p, q). Define 
G(i,r) := 



pk a (t, t) if t < t, 
qk a {r, t) if t > t. 
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For all / G L 2 {[a,b]) one has K%f(t) = j a G{t,T)f(r)dT with G(t,r) G L 2 (A). It 
is not difficult to see that Kp is linear and Kpf G L 2 ([a, b\) for all / G L 2 ([a, b]). 
Moreover, applying the Cauchy-Schwarz inequality and Fubini's theorem, we obtain 



11^/112 = 



< 



G(t,T)f(r)dr 

b 



dt 



\G(t,r)\Urj 
i, ,-b 

\G{t,T)\ 2 drdt. 



\f(r)\ 2 dr 



dt 



For / G L 2 {[a,b]) such that ||/|| 2 < 1 we have \\K%f\\ 2 < (j* j" \G(t,r)\ 2 drdt)' 2 . 

1 

Therefore, \\K^\\ 2 < (f* J* \G(t,r)\ 2 drdt^ 2 . 



□ 



Theorem 2.3. Let k a be a difference kernel, i.e., k a (t,r) = k a (t — r) and k a G 
Li([a,b]). Then Kp : Li([a,b]) — > Li([a,b]) is a well defined bounded linear 
operator. 

Proof. Obviously, the operator is linear. Let a G (0,1), P = (a,t,b,p,q), and 
/ G Li ([a, 6]). Define 

pv ja f P \k a (t-r)\-\f(r)\ if r<t 

1 ' J I q\k a (r-t)\-\f(r)\ ifr>t 

for all (r, t) G A = [a, fe] x [a, b]. Since F is measurable on the square A we have 



/;(/>,*>*) *=/;[i/wi(/ t 



|/(r)| / p|fc a (t-r)|dt + 



J a 



(r - i)| 



< y |/(r)|(p-g)||fe a ||dT 

= (p-«/)IIM- 11/11 • 

It follows from Fubini's theorem that F is integrable on the square A. Moreover, 

11*2/11 = f P f k a (t-T)f(T)dT+q f k a {T-t)f{T)dT dt 
J a J a J t 

< f IP f IM* - r)| • |/(r)| aV + q T |fc Q (r - i)| ■ |/(r)| dr j dt 

J a \ J a J t J 



a \J a 



F{T,t)drj dt 

<(p-«/)IIM -11/11- 

Hence, ifp : L x ([a, 6]) -> Li ([a, 6]) and < (p - q) \\k a \\. 



□ 



Theorem 2.4. Let fci_ a &e a difference kernel, i.e., ki- a (t,r) = ki- a (t — r) and 
ki- a G Li([a, 6]). // / G AC([a,b\), then the K-op of order 1 — a and p- set 
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P = (a,t,b,p,q), i.e., 

t b 
K 1 p - a f(t)=p J h- a (t-T)f(r)dT + q J h- a (T-t)f(T)dr, 

a t 

belongs to AC ([a, b]). 

Proof. Let P 1 = (a,t,b,p,0) and P 2 = (a,t,b,0,q). Then, K x p a = Kp' 01 + Kp' 01 . 
First we show that Kp~ a f G AC ([a, b}). The condition / G AC ([a, b}) implies 

f{x) = j g(t)dt + /(a), where g £ L x ([a, b}) . 

J a 

Let s = x — a and 

h(s) = / k 1 - a (T)g(s + a - r)dr. 
Jo 

Integrating, 

\ h(6)d6= dOl k 1 _ a (T)g(9 + a - r)dr, 
Jo Jo Jo 

and changing the order of integration we obtain 

/>S r-S r-S r-S />S 

/ h(9)d9 = dr k 1 ^ a (r)g{9 + a - T )d9 = I fci_ a (r)rfr / g{6 + a-T)d6. 
Jo Jo Jt Jo Jt 

Putting £ = 9 + a — r and d£ = d9, we have 

P S PS P CC — T 

/ h(9)d9 = / fci_ a (r)dT / 

JO JO Ja 

Because / — /( x — r ) — /( a )i tne following equality holds: 

J a 

f h(6)d6= [ h. a (T)f(x - r)dr - f(a) [ k^ a (r)dr, 

Jo Jo Jo 

that is, 

/ k 1 - a {T)f(x-T)dT= [ h(9)d9 + f(a) f ki- a (T)dT. 

Jo Jo Jo 

Both functions on the right-hand side of the equality belong to AC ([a, b]). Hence, 

f S ki- a (r)f(x-r)dre AC([a,b}). 
Jo 

Substituting t = x — r and dt = —dr, we get 

/ k\- a (x — t)f(t)dt G AC ([a, b\) . 

J a 

This means that Kp^ a f G AC ([a, b}). The proof that Kp~ a f G AC ([a, b}) is 
analogous, and since the sum of two absolutely continuous functions is absolutely 
continuous, it follows that Kp~ a f G AC ([a, b]). □ 
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Remark 2.5. The K-op reduces to the classical left or right Riemann-Liouville 
fractional integral (see, e.g., [25l|39]) for a suitably chosen kernel k a (t,r) and p-set 
P. Indeed, let k a (t - r) = ^(t - r) Q_1 . IfP = (a, i, 6, 1,0), then 

t 

(2.1) ff£/(t) = f L^J(t- TY- l f{r)dT =: a T?f(t) 

a 

is the left Riemann-Liouville fractional integral of order a; if P = (a, t, b, 0, 1), then 

b 

(2.2) Km) = p^y /(r - a_1 /(r)dr =: 

t 

is the right Riemann-Liouville fractional integral of order a. Theorem 12.31 with 
fc a (i — r) = p^y(i — r) a_1 asserts the well-known fact that the Riemann-Liouville 
fractional integrals „if, t l£ : Li ([«,&]) -> £-i ([a, &]) given by (|2.1[) and (|2.2|l arc 
well defined bounded linear operators. 

The fractional derivatives Ap and Bp are defined with the help of the generalized 
fractional integral K-op. 

Definition 2.6 (Generalized Riemann-Liouville fractional derivative). Let P be a 
given parameter set. The operator Ap, < a < 1, is defined by Ap = D o Kp~ a , 
where D denotes the standard derivative. We refer to as operator A (A-op) of 
order a and p-set P. 

A different fractional derivative is obtained by interchanging the order of the 
operators in the composition that defines Ap. 

Definition 2.7 (Generalized Caputo fractional derivative). Let P be a given pa- 
rameter set. The operator Bp, a e (0,1), is defined by Bp = Kp~ a o D and is 
referred as the operator B (B-op) of order a and p-set P. 

Remark 2.8. The operator Bp is defined for absolute continuous functions / G 
AC ([a, b]), while the operator Ap acts on the bigger class of functions / such that 
K p - a fe AC([a,b}). 

Remark 2.9. The standard Riemann-Liouville and Caputo fractional derivatives 
(see, e.g., [251 [3H]) are easily obtained from the generalized operators A p and Bp, 
respectively. Let ki- a (t - r) = r(1 1 _ a) (t - rT", a G (0,1). If P = (a, t, 6, 1,0), 
then 

^/(*) = ^^/(*-r)-/(r)dSr=: aDf /(t) 

a 

is the standard left Riemann-Liouville fractional derivative of order a while 

t 

B a P f(t) = ^^—^ J(t- r)- a f'(r)dT =: ?D«/(t) 

a 

is the standard left Caputo fractional derivative of order a; if P = (a, t, 6, 0, 1), then 

b 

-Am) = - r^s / - *)- a /W* == WO 
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is the standard right Riemann-Liouville fractional derivative of order a while 

b 

t 

is the standard right Caputo fractional derivative of order a. 



3. Main results 

We begin by proving in Section 13.11 that for a certain class of kernels there ex- 
ists a direct relation between the fractional derivatives A P and Bp (Theorem 13. 1|) . 
Section 13.21 gives integration by parts formulas for the generalized fractional set- 
ting (Theorems 13. 4[ 13.71 and l3.8|) . Section 1331 is devoted to variational problems 
with generalized fractional-order operators. New results include necessary opti- 
mality conditions of Euler-Lagrange type for unconstrained (Theorem 13.111) and 
constrained problems (Theorem l3.22p . and a general transversality condition (The- 
orem GTTZ1)- Interesting results are obtained as particular cases. Finally, in Sec- 
tion 13.41 we provide a class of generalized fractional problems of the calculus of 
variations for which one has a coherent embedding, compatible with the least ac- 
tion principle (Theorem 13. 26p . This provides a general answer to an open question 
posed in [13] . 



3.1. A relation between operators A and B. Next theorem gives a useful 
relation between A-op and B-op. In the calculus of variations, equality (I3.2[) can 
be used to provide a necessary optimality condition involving the same operators 
as in the data of the optimization problem (cf. Remark 13.141 of Section 13.31) . 

Theorem 3.1. Let < a < 1, P = (a,t,b,p,q) , and y £ AC ([a, b}). If kernel 
fci_ Q is integrable and there exist functions f and g such that 

(3.1) f kt- a (0,T)d9+ f kt- a (t,0)d0 = g(t) + f(T) 

for all t,r G [a, 6], then the following relation holds: 

(3.2) A P y(t) = py{a)kx- a {t, a) - qy{b)k^ a {b, t) + B^y{t) 
for all t £ [a, b]. 

Proof. Let h\^ a be defined by /ii_ Q (t,r) :— ki- a (t,9)d9 — g(t). Then, by hy- 
pothesis (|3.1[) . c^i-a = —d\hi- a = fci_ Q . We obtain the intended conclusion 
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from the definition of A-op and B-op, integrating by parts, and differentiating: 

,T)y(r)dT + q fci_ a (r, t)y{r)dr \ 



A%y{t) = j^p^yit) = j t jp jf 



-|^< py{t)h\- a {t,T)\ T=t - py(a)h 1 - a (t, a) -p ( h 1 ^ a {t,T)-^y{T)dr 



qy(b)hi- a (b,t) + qy(t)hi- a (T,t)\ T=t + q / hi- a (r, t)— y(r)dr 



d 



py(t)-^hi- a (t,t + e) 



e=0 



qy{b)d 2 h 1 ^ a (b,t) + qy(t)^-hi- a (t + e,t) 



■ py(a)dihi- a (t,a) - p / dihi- a (t,r)—y(r)dT 

b d 
+ q I d2hi- a (T,t)—y(T)dT 



e=0 



= py(o)fcx_„(t,a) - qy(b)ki_ a (b,t) + Bpy(t). 

□ 

Example 3.2. Let k\- a {t — r) = p(i-a) — r )~ Q - Simple calculations show that 
(|3.ip is satisfied. If P = (a, t, b, 1, 0), then (|3.2j) reduces to the relation 

c a D?y{t) = a D?y(t) - ^^—M - a)" Q 

r(i - a) 

between the left Riemann-Liouville fractional derivative Q Z?f and the left Caputo 
fractional derivative „ Df; if P = (a, t, b, 0, 1), then we get the relation 



fD^y(t) = t D^y(t)--A-i 
T(l - a) 



between the right Riemann-Liouville fractional derivative tD^ and the right Caputo 
fractional derivative . 

3.2. Fractional integration by parts. The proof of Theorem 13. II uses one basic 
but important technique of classical integral calculus: integration by parts. In this 
section we obtain several formulas of integration by parts for the generalized frac- 
tional calculus. Our results are particularly useful with respect to applications in 
dynamic optimization (cf. Section [3.3l) . where the derivation of the Euler-Lagrange 
equations uses, as a key step in the proof, integration by parts. 

In our setting, integration by parts changes a given p-set P into its dual P* . The 
term duality comes from the fact that P** = P. 

Definition 3.3 (Dual p-set). Given a p-set P = (a, t, b,p, q) we denote by P* the 
p-set P* — (a, t, b, q,p). We say that P* is the dual of P. 

Our first formula of fractional integration by parts involves the K-op. 

Theorem 3.4 (Fractional integration by parts for the if -op). Let a £ (0,1), 
P = (a,t,b,p,q), k a be a square-integrable function on A = [a, b] x [a, b], and 
f>9 £ L 2 ([a 1 b}). The generalized fractional integral satisfies the integration by 
parts formula 

b b 
(3.3) / g{t)K a P f{t)dt = f f(t)K™,g(t)dt, 
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where P* is the dual of P. 



Proof. Let a G (0,1), P = (a,t,b,p,q), and f,g G L 2 ([a,b}). Define 



r( rt \-l \pUt,r)\-\g(t)\-\f(T)\ ifr<t 
^ lT ' tj -\ \qk a (r,t)\-\g(t)\-\f(r)\ if r > t 



for all (r, t) G A. Then, applying Holder's inequality, we obtain 



| / (r) | ^ \pk a {t,T)\-\g(t)\dt + £\ qk a (r, t)\ ■ \g(t)\ dtj 
f \f(r)\( f b \pka(t,r)\-\g(t)\dt+ f ' \qk a {r,t)\ ■ \g{t)\ dt) 

J a \J a J a I 



< 



dr 



dr 



2 I rb 



\pk a (t,r)\ z dt) / \g(t)\ z dt 



+ ( f \qk a (T,t)\ 2 dt 
l J a 



2 / r b 



\g(t)\ 2 dt 



dr. 



By Fubini's theorem, functions k a . T (t) := k a (t,r) and k a , T (t) := k a (r,t) belong to 
L 2 ([a, b]) for almost all r G [a, b]. Therefore, 



r| l/(T)l (/ |^(*' T )| 2d *) [[Wtfdt 



+ U\qk a {r,t)\ 2 dt) (J*\9(t)\ 



dt 



dr 



= IMI 2 j[ [l/Wl(l|P*a,rll2+||9*a,r|| 2 ) 

<\\9h^j\f(r)\ 2 dT^ ^ 6 |||pfc a;T || 2 + 
< IMI 2 ' H.fll 2 (lbM 2 + Il9*a|| 2 ) < 00. 



dr 



qk 



dr 



GENERALIZED FRACTIONAL CALCULUS WITH APPLICATIONS 



9 



Hence, we can use again Fubini's theorem to change the order of integration: 

b b t b b 

g{t)K£f(t)dt=p Jg(t)dt J f(T)k a (t,r)dT + q J g(t)dt J /(r)fc Q (r, t)dr 

a a at 

6 6 6 r 

p J f(r)dT J g{t)k a {t,T)dt + q J f(r)dr J g(t)k a ( T ,t)dt 

a t a a 

b 

f(r)K^g(r)dT. 

□ 

Next example shows that one cannot relax the hypotheses of Theorem 13.41 

Example 3.5. Let P = (0,t, 1, 1, -1), f{t) = g{t) = 1, and k a {t,r) = ( fcf )a . 
Direct calculations show that 

* t 2 -r 2 . r 1 r 2 -t 2 



dr ) dt = / dt = - 



and 



o (t 2 +r 2 ) 2 J J t 2 + l 4 



T r 2 -t 2 , f 1 t 2 -r 2 



*V*r= M-/ o w^f dt + j T WT^Y 2dt ) dT 

dt \ dr = — / — = dr 



/o (t 2 + r 2 ) 2 ; y r 2 + i 4' 

Therefore, the integration by parts formula (|3.3[) does not hold. Observe that in 
this case L \k a (t, r)| 2 drdt — oo. 

For the classical Riemann-Liouville fractional integrals the following result holds. 

Corollary 3.6. Let \ < a < 1. If f,g G L 2 ([a,b]), then 

(3.4) f g{t) J?f{t)dt = f ' f(t) t I«g(t)dt. 

J a J a 

Proof. Let k a (t,r) = p^y(^ — r) Q_1 . For a G (5, l) , fc a is a square-integrable 
function on A (see, e.g., (221 Theorem 4]). Therefore, (|3.4j) follows from (|3.3|) . □ 

Theorem 3.7. Let < a < 1 and P = (a,t,b,p, q) . If k a (t,r) — k a (t — r), 
k a G L\{[a,b\), and f,g G C ([a, b]), then the integration by parts formula 
holds. 

Proof. Let a G (0,1), P — (a,t,b,p, q) , and /,g G C([a,b]). Define 

MrA-t ^M*-T)|-Ifl(*)|-|/Wl 
^ T '*J- 1 |g*a(r-*)|.|ff(*)|-|/W| ifr>* 
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for all (t, t) E A. Since / and g are continuous functions on [a, b], they are bounded 
on [a, b], i.e., there exist real numbers Ci, C2 > such that |<?(t)| < C\ and |/(£)| < 
C2 for all i E [a, 6]. Therefore, 



^ ^ F{r,t)dtj dr 



|/(r)| / |p*„(*-r)|. ]«/(*)! dt+ / |g*a(r-*)|-b(*)|d« 



< 



l/MI / |pfca(*-T)|-|5(*)l*+ / \qk a (r-t)\-\g(t)\dt 



dr 



dr 



< C1C2 



C1C2 



a \ ^ a 



\pk a (t - t)\ dt + \qk a (r -t)\dt\ dr 



\q\) \\k a \\ (b-a) < 00. 



Hence, we can use Fubini's theorem to change the order of integration in the iterated 
integrals. □ 

The next theorem follows from the classical formula of integration by parts and 
fractional integration by parts for the K-op. 

Theorem 3.8 (Fractional integration by parts for A-op and B-op). Let a E (0, 1) 

and P — (a, t, b,p, q). If f,g E AC {[a, b]), then 

b 



(3.5) 



g(t)A a P f(t)dt = g(t)K' p - a f(t) 



f(t)B$.g{t)dt, 



(3.6) 



g{t)B a P f{t)dt = f(t)K p -, a g(t) 



f(t)A a P ,g(t)dt. 



Proof. From Definition UM one has A P f(t) = DKf a f{t). Therefore, 



g(t)A a P f(t)dt= I g(t)DK p - a f(t)dt 



g(t)K p - a f(t)\ 



Dg(t)K p - a f(t)dt, 



where the second equality follows by the standard integration by parts formula. 
From (13. 3[) of Theorem 13.41 it follows the desired equality (|3 . 5|) : 



g(t)A a P f(t)dt= g(t)K p - a f(t)\ b a - / f(t)K p - a Dg(t)dt 



We now prove flM]). From Definition O we know that B£f(t) = Kp- a Df{t). It 
follows that 



By Theorem 13.4 



f g(t)B$f(t)dt = f g(t)K p - a Df(t)dt. 
f g{t)B a P f{t)dt = f Df(t)K p - a g(t)dt. 

J a J a 
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The standard integration by parts formula implies relation (|3.6|) : 

' g(t)By(t)dt= f(t)K 1 p - a g(t)\ b a - [ f(t)DK p -, a g(t)dt 



Corollary 3.9. Let < a < 1. If f,g £ AC([a,b}), then 

g{t) c a D?f{t)dt = f(t)jl- a g(t)t b + I f(t) t D%g(t)dt 



□ 



t—a 

b , r b 

t=b 



f(t) a D?g(t)dt= f(t) a I t i - a g(t)\ t=a + / g{t)\DUit)dt. 

J a 

3.3. Fractional variational problems. We study variational functionals with 
a Lagrangian depending on generalized Caputo fractional derivatives as well as 
derivatives of integer order. Note that the only possibility of obtaining y' from 
Bpy or Apy is to take the limit when a tends to one but, in general, such a limit 
does not exist [U]. Moreover, our Lagrangians may also depend on generalized 
fractional integrals. This last possibility is used in Section [3T4l to solve the important 
coherence problem. 

Our proofs are easily adapted to the cases when one considers Riemann-Liouville 
A-op derivatives instead of Caputo _B-op derivatives, and vector admissible func- 
tions y instead of scalar ones. Such versions are left to the reader. 

3.3.1. Fundamental problem. We consider the problem of extremizing (minimizing 
or maximizing) the functional 

b 

(3.7) J[y} = J F(t,y(t),y'(t),B^y(t),K P2 y(t))dt 

a 

subject to boundary conditions 

(3-8) y(a) = y a , y(b) = y b , 

where a, /3 £ (0, 1) and Pj — (a, t, b,pj,qj), j = 1, 2. 

Definition 3.10. A Lipschitz function y £ Lip ([a, b]; R) is said to be admissible 
for the fractional variational problem (|3.7[) - (|3.8[) if it satisfies the given boundary 
conditions (13.81). 



For simplicity of notation we introduce the operator {■}pfp 2 defined by 
Mpfp 2 (*) = (t,y(t),y'(t),B^y(t),K^y(t)" 



We can then write (|3.7p in the form J[y] = J F {y}p'fp 2 (t)dt. We assume that 

a 

F £ C 1 ([a, b] x M 4 ; R) , t M> d^F {y} c pfp 2 (t) is absolutely continuous and has a 
continuous derivative A P » , and t H> d$F {y}pfp 2 (t) has a continuous derivative 

d_ 

dt " 

Next result gives a necessary optimality condition of Euler-Lagrange type for 
problem (pTTf - tprgl) . 
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Theorem 3.11. Let y be a solution to problem (|3.7|) - (|3.8j) . Then, y satisfies the 
generalized Euler-Lagrange equation 

(3.9) d 2 F{y} a p f P2 (t) - ±d 3 F{y} a p f P2 (t) - A P ,d,F {y} a p f P2 (t) 

+ K^d 5 F{y} a p f P2 (t)=0 

for t e [a, b] . 

Proof. Suppose that y is an extremizer of J . Consider the value of J at a nearby 
admissible function y(t) = y(t) + erj(t), where e G K is a small parameter and 
r\ G Lip ([a, b]; K) is an arbitrary function satisfying 77(a) = rj(b) = 0. Let J(e) :— 

b 

S\y] = f F {y + £Tj}'pfp 2 (t)dt. A necessary condition for y to be an extremizer is 

a 

given by J'(O) = 0, i.e., 

(3.10) J id 2 F{y} p f P2 (t) ■ V (t) + d 3 F{y} a p f P2 (t)-rj(t) 

a ^ 

+ d,F {y} a p f P2 (t)BpMt) + d 5 F {y} a p f P2 (t) ■ K^{t) J dt = 0. 

Using the classical integration by parts formula as well as our generalized fractional 
versions (Theorems 13.41 13.71 and I3.8j) we obtain that 

6 b 

J diFB^ndt = - J r]A P ,d 4 Fdt + r]Kp- ; a d A F 

a a 

and 

t 

d 5 FK p rjdt = I r]K p ,d 5 Fdt, 



where P* = {a,t 7 b,qj,pj), j = 1,2, is the dual of Pj. Because -q(a) = r](b) = 0, 
(|3.10|) simplifies to 

£ V (t) (d 2 F {y} a p f P2 (t) - j t d,F {y} p f P2 (t) A P ,d,F {y} a p f P2 (t) 

+ K Pi d 5 F{ y } a p f P2 (t)\dt = 0. 

We obtain (|3.9|) applying the fundamental lemma of the calculus of variations. □ 

Remark 3.12. If the functional (|3.7[) does not depend on Kp 2 y{t) and B Pi y(t), then 
Theorem 13.111 reduces to the classical result: if y is a solution to the problem 
6 

F (t, y(t), y'(t)) dt — > extr, y(a) = y a , y(b) = y b) 



GENERALIZED FRACTIONAL CALCULUS WITH APPLICATIONS 



13 



then y satisfies the Euler-Lagrange equation 

(3.11) d 2 F (t, y{t), y'{t)) - j^F (t, y(t),y'(t)) = 0, 
t G [a, b]. 

Remark 3.13. In the particular case when functional (|3.7p does not depend on the 
integer derivative of function y, we obtain from Theorem 13.111 the following result: 
if y is a solution to the problem of extremizing 

b 

J[y] = J F (t,y{t),B^y{t),K^y{tj) dt 

a 

subject to y(a) = y a and y(b) — yb, where a,ft£ (0,1) and Pj — (a,t,b,Pj,qj), 
j = 1,2, then d 2 F - A P ,d 3 F + K^,d 4 F = with P* = (a,t,b,qj,pj), j = 1,2. 
This extends some of the recent results of [4] . 

Remark 3.14. The optimization problem (|3.7[) - (I3.8[) does not involve the general- 
ized Riemann-Liouville fractional derivative A-op while the necessary optimality 
condition (I3.9|) does. However, using Theorem 13 . 1 1 the Euler-Lagrange equation 
(|3.9[) can be written in terms of B-op as 

(*) - |Wy}^ 2 (t)-qd,F{y} a p f P2 (a)h^ a (t,a) 
+ pd 4 F {y} a p f P2 (b) t) - Bp* d±F {y} P f Pj (t) + K P ,d 5 F {y} a p f P2 (t) = 0. 

Corollary 3.15. Let < a, ft < 1. If y is a solution to the problem 

b 

[ F (t, y{t), y'(t), C a D«y{t), J t P y(t)) dt — > min , y(a) = y a , y(b) = y b , 

J \ / y£Lip 

a 

then the following Euler-Lagrange equation holds: 

t D<Zd A F (t, y(t),y'(t)C D?y(t), a f t y{t)) + t f b d 5 F (t, y(t),y'(t)° Dfy(t), a f t y{t)) 

+d 2 F (t, y(t),y'(t)C D?y(t), a f t y{t))-jd^F (t, y(t),y'(t)^ D?y(t), a f t y[t)) = 0. 

Proof. The intended Euler-Lagrange equation follows from (|3.9p by choosing the 
p-sets P\= P 2 = (a, t 7 b, 1, 0) and the kernel ki_ a (t — r) = r (i_ a ) (t — r)~ a . □ 

Corollary 3.16. Let J be the functional 
b 

J[y] = J F(t,y(t),y'(t),P^Dfy(t) + qfD^y(t))dt, 

a 

where p and q are real numbers, and y be an extremizer of J satisfying boundary 
conditions y(a) = y a and y{b) = yb- Then, y satisfies the Euler-Lagrange equation 

(3.12) p tDfdiF + q a D?d 4 F + d 2 F - ^-d 3 F = 

dt 

with functions evaluated at (t, y(t), y'(t),p^D^y(t) + qf D^y(t)^j , t £ [a, b\. 
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Proof. Choose Pi = (a,t,b,p, —q) and kx- a (t — r) = r(i-a) (t — T )~ a ■ Then the 
B-op reduces to the sum of the left and right Caputo fractional derivatives and 
(l3~T21) follows from flJSJ}. □ 



3.3.2. Free initial boundary. Let in problem (|3.7|) - (|3.8|) the value of the unknown 
function y be not preassigned at the initial point t = a, i.e., 

(3.13) y(a) is free and y(b) = yb- 

Then, we do not require 77 in the proof of Theorem 13. Ill to vanish at t — a. There- 
fore, following the proof of Theorem 13. Ill we obtain 

(3.14) n(a)d 3 F {y} P f P2 (a) + 77(a) K% a d 4 F {y}£? ft (t) 

1 t—a 

+ j\(t){d 2 F{y} a Pi P P2 (t) - ±d 3 F{y} p f P2 (t) - A P ,d,F {y} a p f P2 (t) 

+ K P [,d 5 F{y} a p f P2 (t)y t = 

for every admissible 77 £ Lip([a, b];M.) with 77(6) = 0. In particular, condition (|3.14j) 
holds for those n that fulfill 77(a) = 0. Hence, by the fundamental lemma of the 
calculus of variations, equation (|3.9p is satisfied. Now, let us return to (13. 14)) and let 
77 again be arbitrary at point t = a. Using equation (|3.9p . we obtain the following 
natural boundary condition: 



(3.15) d 3 F {y} a p f P2 (a) + K p ^d 4 F {y}^ (t) 

We just obtained the following result. 



= 0. 

t — a 



Theorem 3.17. If y 6 Lip([a, b];M.) is an extremizer of (|3.7j) subject to y(b) = yb, 
then y satisfies the Euler-Lagrange equation (|3.9p and the transversality condition 

Corollary 3.18 (cf. Theorem 2.3 of [9]). Let J be the functional given by 

b 

J[y] = f F(t,y(t)CD?y(t))dt. 



If y is a local minimizer of J satisfying the boundary condition y(b) — yb, then y 
satisfies the Euler-Lagrange equation 

(3.16) W (t,y®£ D?v(t)) +t ^d 3 F (t,y(t)° D?y(t)) = 

and the natural boundary condition 

(3-17) t ll- a d 3 F(t,y(t)^D?y(t))\ t=a = 0. 



Proof. Let functional Q3.7P be such that it does not depend on the integer derivative 
y'(t) and on K-op. If Pi = (a, t, b, 1, 0) and ki- a (t — r) = m— a) ft — r )~ Q ; then 
B-op reduces to the left Caputo fractional derivative and from (|3.9p and (|3.15[) we 
deduce (|3.16p and (|3.17p . respectively. □ 
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Remark 3.19. Observe that if the functional (|3.7[) is independent of K-op, then the 
problem defined by (|3.7|) and (|3.13[) takes the form 

b 

J F (t, y(t),y'(t), Bp ± y(t)) dt — >• extr, y{b) = y b 

a 

(y(a) free) and the optimality conditions (13.91) and (I3.15P reduce respectively to 

d 2 F (t, y(t), y'(t), B a Pi y{t)) - j^F (t, y(t),y'(t),B^y(t)) 

-A P .d 4 F(t,y(t),y'(t),B p - 1 y(t))=0 



SDdd 3 F{a,y{a),y'{a),B% L y(a)) + K^ a d A F (t,y(t),y'(t),B^y(t)) 



0. 



3.3.3. Isoperimetric problems. One of the earliest problems in geometry was the 
isoperimetric problem, already considered by the ancient Greeks. It consists to 
find, among all closed curves of a given length, the one which encloses the maxi- 
mum area. The general problem for which one integral is to be given a fixed value, 
while another is to be made a maximum or a minimum, is nowadays part of the 
calculus of variations [341 135) . Such isoperimetric problems have found a broad 
class of important applications throughout the centuries, with numerous useful im- 
plications in astronomy, geometry, algebra, analysis, and engineering [TUJITB]. For 
recent advancements on the study of isoperimetric problems see [SJ [T7] and refer- 
ences therein. Here we consider isoperimetric problems with generalized fractional 
operators. Similarly to Sections 13.3.11 and 13.3. 2\ we deal with integrands involving 
both generalized Caputo fractional derivatives and generalized fractional integrals, 
as well as the classical derivative. 

Let < a, j3 < 1 and Pj = (a, t, b,Pj, qj), j — 1, 2, be given p-sets. Consider the 
following isoperimetric problem: 

6 

(3.18) J[y] = J F {y} P f P2 (t)dt — > extr, 

a 

(3-19) y(a)=y a , y(b) = y b , 

b 

(3.20) X[y} = J G{y} P fp 2 (t)dt = t 

a 

Definition 3.20. A Lipschitz function y : [a, b] — > K is said to be admissible for 
problem (|3.18|l - (|3.20|) if it satisfies the given boundary conditions (|3.19|) and the 
isoperimetric constraint (|3.20[) . 

We assume that F,G € C 1 ([a, b] x M , R), £ is a specified real constant, func- 
tions t H> dt±F {y}pf p 2 (t) and t H> d^G {y} p ' p (t) are absolutely continuous 
and have continuous derivatives Ap, , and functions t h-> <9 3 F {y}pfp 2 (t) and 
t i — y {y} P fp 2 {t) have continuous derivatives ^. 
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Definition 3.21. An admissible function y € Lip([a,b],M.) is an extremal for 1 if 
it satisfies the Euler-Lagrange equation p.9[) associated with (|3.20|) . i.e., 

d 2 G {y} a p f P2 (t) - ±d 3 G {y} a p f P2 (t) A P ,d,G {y} a p f P2 (t) 

where P* = (a,t,b,qj,Pj), j = 1,2, and t G [a,b]. 

The next theorem gives a necessary optimality condition for the generalized 
fractional isoperimetric problem (|3.18p - (|3.20l) . 

Theorem 3.22. If y is a solution to the isoperimetric problem (|3.18[) - (|3.20l) and 

is not an extremal for I, then there exists a real constant A such that 

(3.21) d 2 H{y} a p f P2 (t) - ±d 3 H{y} p f P2 (t) A%.d*H {y}^ (t) 

+ K%d 5 H{y} p f P2 (t) = 0, 

t G [a, b], where H{t, y, u, v, w) = F(t, y, u, v, w) — XG(t, y, u, v, w). 

Proof. Consider a two-parameter family of the form y = y + £if]i + £ 2 r\ 2 , where for 
each i G {1, 2} we have %(a) = r)%(b) = 0. First we show that we can select rj 2 such 
that y satisfies (|3.20l) . Consider the quantity 

6 

T[y] = j gU y(t) + e im {t) + e 2 n 2 {t), j t (v(*) + + 9»»ft(t)) , 

a 

B Pl (y(t) + £iVi{t) + £2V2{t)) , K P P2 (y(t) + s lVl (t) + £ 2 V2{t))jdt. 

Looking to X[y] as a function of £\,£2, we define j(ei,£ 2 ) : = 1[y] — Thus, 
1(0, 0) = 0. On the other hand, integrating by parts, we obtain 



dl 

de 2 



V2(t)(d 2 G{y} a p f P2 (t) - j t d 3 G{y} a P f jP2 (t) 



(0,0) 

- A P ,d A G {y} a p f P2 (t) + K p *d 5 G {y}^ (t) ) dt, 



where P* = (a,t,b,qj,pj), j — 1,2. We assumed that y is not an extremal for I. 
Hence, the fundamental lemma of the calculus of variations implies that there exists 

a function ri 2 such that M^- 7^ 0. According to the implicit function theorem, 

062 (0,0) 

there exists a function £2(0 defined in a neighborhood of such that I(£\, £ 2 (£i)) = 
0. Let J(ei, £2) = J[y\- Function J has an extremum at (0, 0) subject to 1(0, 0) = 0, 
and we have proved that V/(0, 0) ^ 0. The Lagrange multiplier rule asserts that 
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there exists a real number A such that V(J(0,0) — AI(0, 0)) = 0. Because 
b 



dJ 
dei 



(0,0) 



I (d 2 F{y} a p f P2 (t) - ±d 3 F{ V r p f P2 (t) 

a 

- A P » d,F {y} a p f P2 (t) + KL d 5 F {y} a p f P2 (t) ) m (t)dt 



and 



dl 

dei 



(0,0) 



(9 2 G{yr p f P2 (t)-j t d 3 G{y} P f tP2 (t) 

- A P ,d 4 G{y} a P f P2 (t) + K^d 5 G{y} a P f P2 (t) \ Vl (t)dt, 



one has 

6 



J [d 2 H{y} p f P2 (t) - ±d 3 H{y} p f P2 (t) 



A P ,d A H{ y } a p f P2 (t)+K p ,d 5 H{y} p ^ P2 (t) \ Vl (t)dt = 



a,/3 



We get equation (|3 . 2 1 1) from the fundamental lemma of the calculus of variations. 

□ 



As particular cases of our problem (|3.18[) - (13.20[) . one obtains previously studied 
fractional isoperimetric problems with Caputo derivatives. 

Corollary 3.23 (cf. Theorem 3.3 of Let y be a local minimizer to 
J[y] = L(t, y(t), C a D?y{t)) dt — ► min, 



l[y]= / g(t,y(t) > %D<Xy(t))dt = Z, 

y(a) = y a , y(b) = yb- 

If y is not an extremal of I, then there exists a constant A such that y satisfies 
d 2 F(t,y{t), c a Dfy{t))+ t D^d 3 F(t,y{t), c a Dfy{t))^Q,te[a,b], withF = L+\g. 



Proof. The result follows from Theorem 13. 221 by choosing the kernel ki^ a (t — r) = 
r(i-a) — 7") _Q and the p-set Pi to be Pi — (a,t,b, 1,0). Indeed, in this case the 
operator —A P , becomes the right Riemann-Liouville fractional derivative, and the 
operator B P becomes the left Caputo fractional derivative. □ 



Remark 3.24. If functionals (|3.18j) and (|3.20j) do not depend on integer derivatives, 
then problem (|3.18p - (|3.20[) is reduced to extremize functional 



J[y) = J F {t,y(t),B^y(t),K^y(t)) dt 
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subject to boundary conditions y(a) = y a , y{b) — yb, and the isoperimetric con- 
straint 

b 

X[y] = J G(t,y(t),B^y(t),K p P2 y(t))dt. 

a 

By (|3~2Tj) there exists A such that y satisfies d 2 H-A P ,d i H+K P '»d 5 H = Q,te [a, b], 
with H = F-XG. 

Remark 3.25. Theorem 13.221 can be extended to the case when y is an extremal for 
X. The proof is similar but one needs to use the extended (abnormal) Lagrange 
multiplier rule. The method is given in [6]. 

3.4. The coherence embedding problem. The notion of embedding introduced 
in [14] is an algebraic procedure providing an extension of classical differential equa- 
tions over an arbitrary vector space. This formalism is developed in the framework 
of stochastic processes [14] , non-differentiable functions [15] , discrete sets [H] , and 
fractional equations [13]. The general scheme of embedding theories is the fol- 
lowing: (i) fix a vector space V and a mapping i : C°([a,b],M. n ) — > V; (ii) extend 
differential operators over V; (iii) extend the notion of integral over V. Let (i, D, J) 
be a given embedding formalism, where a linear operator D : V — > V takes place 
for a generalized derivative on V, and a linear operator J : V — > R takes place for 
a generalized integral on V. The embedding procedure gives two different ways, a 
priori, to generalize Euler-Lagrange equations. The first (pure algebraic) way is to 
make a direct embedding of the Euler-Lagrange equation. The second (analytic) 
is to embed the Lagrangian functional associated to the equation and to derive, by 
the associated calculus of variations, the Euler-Lagrange equation for the embed- 
ded functional. A natural question is then the problem of coherence between these 
two extensions: 

Coherence problem. Let [t, D, J) be a given embedding formalism. Do we 
have equivalence between the Euler-Lagrange equation which gives the direct em- 
bedding and the one received from the embedded Lagrangian system? 

For the standard fractional differential calculus of Riemann-Liouville or Caputo, 
the answer to the question above is known to be negative. For a gentle explanation 
of the fractional embedding and its importance, we refer the reader to [321 l23j 
129] . Here we propose a coherent embedding in the framework of our fractional 
generalized calculus by choosing the generalized fractional operator to be Kp with 
q = —p. A direct embedding of the classical Euler-Lagrange equation (|3 . 1 1 [) gives 

(3.22) d 2 F (i, y(t), K$y{t)) - K£d 3 F (t, y(t), K%y{t)) = 

for t E [a,b]. On the other hand, we can apply Theorem 13.111 to the embedded 

b 

Lagrangian functional J[y] = J F (t,y(t), Kpy(t)) dt. Let P = (a,t,b,p, —p) and 

a 

a G (0, 1). If y is a solution to the problem 

b 

(3 23) J[y] = J F{t, y(t), K%y{t)) dt — > extr, 

a 

y(a) = y a , y{b) = yb, 
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then, by Theorem 13.111 y satisfies the Euler-Lagrange equation given by 

(3.24) d 2 F (t, y(t), K%y(t)) + K£,d 3 F (t, y(t),K£y(t)) = 0, 

t G [a, b]. For an arbitrary kernel k a , an easy computation shows that for p = —q 
one has Kpf(t) = —Kp,f(t). Therefore, equation (I3.24[) can be written in the 
form 

(3.25) 3 2 F (t, y(t),K^y(t)) - Kpd^F (t, y(t), K%y(t)) = 0, 

t G [a, b]. It means that the Euler-Lagrange equation (|3.22|) obtained by the direct 
fractional embedding procedure and the Euler-Lagrange equation (|3.25|) obtained 
by the least action principle coincide. We just proved the following result. 

Theorem 3.26. Let k a (t,r) be an arbitrary kernel and P a p-set with q = —p: 
P = (a, t, b,p, —p). Then the fractional variational problem (|3.23l) is coherent. 

4. Illustrative examples 

In this section we illustrate our results through two examples of isoperimetric 
problems with different kernels. Explicit expressions for the minimizers are given. 

In Example 14.11 we make use of the Mittag-Leffler function of two parameters. 
Let a, f3 > 0. We recall that the Mittag-Leffler function is defined by 

\. k 

This function appears naturally in the solution of fractional differential equations, as 
a generalization of the exponential function |12j . Indeed, while a linear second order 
ordinary differential equation with constant coefficients presents an exponential 
function as solution, in the fractional case the Mittag-Leffler functions emerge [25] . 

Example 4.1. Let a G (0, 1) and Consider the following problem: 

J(y) = [ (y' + TTpyf dt — > min, 
Jo 

(4.1) X(y)= f {y'+B P y)dt = ^ 

Jo 

1/(0) = 0, y(l) = J E^ aA (- (1 - t) 1 -") idr, 

where k\^ a (t — r) = r(i-a) ft ~ "r)~ a and P = (0, t, 1, 1, 0). In this case the _B-op 
becomes the left Caputo fractional derivative, and the augmented Lagrangian H of 
Theorem l3.22l is given by H(t, y, v, w) = (v + w) 2 — X(v + w). One can easily check 
that 

(4.2) y(t) = J q Ei-ci (- (t - r) 1 -") £dr 

is not an extremal for I and satisfies y' + B P y = £. Moreover, (|4.2j) satisfies ()3.2ip 
for A = 2^, i.e., 

~ (2 („' + B P y) - 20 - A%. (2 (y' + B P y) - 2f) = 0, 

where P* = (0,t, 1,0, 1) is the dual p-set of P. We conclude that (|4.2j) is an 
extremal for problem (|4.1|) . Since in this example one has a problem (|3.18j) - ()3.20|) 
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with F(t, y, v, w) = (v + w) 2 and G(t, y,v,w) = v + w, simple convexity arguments 
show (see [7j Section 6] and Section 3.4]) that (|4.2p is indeed the global minimizcr 
to problem (|4.1[) . 

Example 4.2. Let a 6 (0, 1), P = (0,t, 1, 1,0). Consider the following problem: 
J(y) = f {K a P y + tf dt — > min, 



/ tK a P ydt = Z, 
Jo 



2/(0) = £ - 1, 2/(1) = (£ - 1) I 1 + J r a (l - r)dr 

where the kernel is such that k a (t, r) = k a (t ~ r) with fc Q (0) = 1 and K P „i 

0. The resolvent r a {t) is given by r a (i) = ^p^y — 1 , k a (s) — C[k a (t)], 

where C and £ _1 are the direct and inverse Laplace transforms, respectively. Since 
Kp,t ^ 0, there is no solution to the Euler-Lagrange equation for functional X. The 
augmented Lagrangian H of Theorem l3.22l is given by H(t, y, w) = (w + t) 2 — Xtw. 
Function 

t 



y(i) = (£-l) (l + jf r Q (i-r)dr) 



is the solution to the Volterra integral equation of the first kind Kpy = (£ — l)t 
(see, e.g., Eq. 16, p. 114 of [10]) and for A = 2£ satisfies our optimality condition 

ma-. 

(4.3) K£, (2 (Kpy + 1) - 2£t) = 0. 

The solution of (14.31) subject to the given boundary conditions depends on the 
particular choice for the kernel. For example, let k a (t — r) = e a ^ T \ Then the 
solution of (|4.3j) subject to the boundary conditions y(0) = £ — 1 and j/(l) = 
(£ - 1)(1 - a) is y(t) = (f - 1)(1 - at) (cf. 40 s p. 15]). If fc a (t - r) = cos (a(t - r)), 
then the boundary conditions are y(0) — £ — 1 and 2/(1) = (£ — 1) (l + a 2 /2), and 
the extremal is y(t) = (£ - 1) (l + a 2 t 2 / 2 ) (cf- 03 P- 46]). 

Borrowing different kernels from book [40j . many other examples of dynamic 
optimization problems can be explicitly solved by application of the results of Sec- 
tion [3j 
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